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Abstract. Improving our understanding of the earliest stages of star 
formation is crucial to gain insight into the origin of stellar masses, multi- 
ple systems, and protoplanetary disks. We discuss recent advances made 
in this area thanks to detailed mapping observations at infrared and 
(sub)millimeter wavelengths. Although ambipolar diffusion appears to 
be too slow to play a direct role in the formation of dense cores, there is 
nevertheless good evidence that the gravitational collapse of isolated pro- 
tostellar cores is strongly magnetically controlled. We also argue that the 
beginning of protostellar collapse is much more violent in cluster-forming 
clouds than in regions of distributed star formation. 
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1. INTRODUCTION 

One of the main limitations in our present understanding of 
the star formation process is that we do not know well the initial 
conditions for protostellar collapse. In particular, there is a ma- 
jor controversy at the moment between two schools of thought for 
the formation and evolution of dense cores within molecular clouds: 
The classical picture based on magnetic support and ambipolar dif- 
fusion (e.g. Shu et al. 1987, 2003; Mouschovias & Ciolek 1999) 
has been seriously challenged by a new, more dynamical picture, 



2 



Ph. Andre, A. Belloche, P. Hennebelle, D. Ward-Thompson 



which emphasizes the role of supersonic turbulence in supporting 
clouds on large scales and generating density fluctuations on small 
scales (e.g. Klessen et al. 2000; Padoan & Nordlund 2002). Improv- 
ing our knowledge of the earliest stages of star formation is thus of 
prime importance, especially since there is now good evidence that 
these stages control the origin of the stellar initial mass function (e.g. 
Motte, Andre, & Neri 1998; Testi & Sargent 1998; Johnstone et al. 

2000) . 

Here, we discuss recent observational advances concerning the 
density and velocity structure of cloud cores (§2 and §3), and we 
compare the results with the predictions of theoretical models (§4). 

2. DENSITY STRUCTURE OF PRESTELLAR CORES 

Two techniques have been used to trace the density structure of 
cloud cores (see Fig. la): (1) mapping the optically thin (sub) millimeter 
continuum emission from the cold dust contained in the cores, and (2) 
mapping the same cold core dust in absorption against the background 
infrared emission (originating from warm cloud dust or remote stars) . 

Ward-Thompson et al. (1994, 1999) and Andre et al. (1996) em- 
ployed the first approach to probe the structure of prestellar cores 
(see also Shirley et al. 2000). Under the simplifying assumption of 
spatially uniform dust temperature and emissivity properties, they 
concluded that the radial density profiles of isolated prestellar cores 
were not consistent with the single p(r) oc r~ 2 power law of the singu- 
lar isothermal sphere (SIS) but were flatter than p(r) oc r~ l in their 
inner regions (for r < i?// at ), and approached p(r) oc r~ 2 only beyond 
a typical radius R f i at ~ 2500-5000 AU. 

More recently, the use of the absorption approach, both in the 
mid-IR from space (e.g. Bacmann et al. 2000, Siebenmorgen & 
Kriigel 2000) and in the near-IR from the ground (e.g. Alves et al. 

2001) , made it possible to confirm and extend the (sub) millimeter 
emission results, essentially independently of any assumption about 
the dust temperature distribution. 

The typical column density profile found by these emission and 
absorption studies of prestellar cores has the following characteristics 
(see, e.g., Fig. lb): 

a) a flat inner region (of radius R f i at = 5000 ± 1000 AU for L1689B 
according to the fitting analysis shown in Fig. lb) , b) a region roughly 
consistent with N n2 oc f" 1 (corresponding to p oc r~ 2 for a spheroidal 
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Fig. 1. (a) ISOCAM 6.75 p,m absorption image of the prestellar core 
L1689B in Ophiuchus (scale on the right in MJy/sr). The 1.3mm con- 
tinuum emission map of Andr et al. (1996) with the IRAM 30 m telescope 
is superposed as contours (levels: 10, 30, 50mJy/13"-beam). (b) Column 
density profile of L1689B (crosses) derived from the absorption map shown 
in (a) by averaging the intensity over elliptical annuli for a 40° sector in 
the southern part of the core. The dashed curves show the most extreme 
profiles compatible with the data given the uncertainties affecting the ab- 
sorption analysis. The solid curve is the best fit of a Bonnor-Ebert sphere 
model (embedded in a medium of uniform column density), obtained with 
the following parameters: pd Pout = 40 ± 15 (i.e., well into the unstable 
regime), T e ff = 50 ± 20 K, P ex tl^B = 5 ± 3 x 10 5 K cm -3 . For comparison, 
the dotted line shows the iV H2 oc f~ l profile of a SIS at T = 10 K. (Adapted 
from Bacmann et al. 2000.) 

core), c) a sharp edge where the column density falls off more rapidly 
than iV H2 oc f -2 with projected radius defining the core outer radius 
(R mt = 28000 ± 1000 AU for L1689B). 

2.1. Comparison with models of core structure 

The results summarized above set strong constraints on the den- 
sity structure at the onset of protostellar collapse. The circularly- 
averaged column density profiles can often be fitted remarkably well 
with models of pressure-bounded Bonnor-Ebert spheres, as first demon- 
strated by Alves et al. (2001) for B68. This is also the case of L1689B 
as illustrated in Fig. lb. The quality of such fits shows that Bonnor- 
Ebert spheroids provide a good, first order model for the structure of 
isolated prestellar cores. In detail, however, there are several prob- 
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lems with these Bonnor-Ebert models. First, the inferred density 
contrasts (from center to edge) are generally larger (i.e., £ 20-80 - 
see Fig. lb) than the maximum contrast of ~ 14 for stable Bonnor- 
Ebert spheres. Second, the effective core temperature needed in 
these models (for thermal pressure gradients to balance self-gravity) 
is often significantly larger than both the average dust temperature 
measured with ISOPHOT (e.g. Ward-Thompson et al. 2002) and the 
gas temperature measured in NH3 (e.g. Lai et al. 2003) . In the case 
of L1689B, for instance, the effective temperature of the Bonnor- 
Ebert fit shown in Fig. lb is T e ff ~ 50 K, while the dust tempera- 
ture observed with ISOPHOT is only T d ~ 11 K (Ward-Thompson 
et al. 2002). Third, the physical process responsible for bounding 
the cores at some external pressure is unclear. These arguments sug- 
gest that prestellar cores cannot simply be described as isothermal 
hydrostatic structures and are either already contracting (e.g. Lee, 
Myers, Tafalla 2001) or experiencing extra support from static or 
turbulent magnetic fields (e.g. Curry & McKee 2000). It has also 
recently been pointed out that good Bonnor-Ebert fits can often be 
found for non-equilibrium, transient "cores" produced by turbulence 
(Ballesteros-Paredes et al. 2003). 

As shown by Bacmann et al. (2000), one way to account for large 
density contrasts and high effective temperatures is to consider mod- 
els of cores initially supported by a static magnetic field and evolving 
through ambipolar diffusion (e.g. Ciolek & Mouschovias 1994, Basu 
& Mouschovias 1995). A serious problem, however, with such models 
is that they need a strong field in the low-density ambient cloud (~30- 
100 fiG) to explain the observed sharp edges (see Bacmann et al. 
2000). This seems inconsistent with existing Zeeman measurements 
(e.g. Crutcher 1999, Crutcher & Troland 2000). Furthermore, when 
the ambient field is strong (highly subcritical case) , the model cores 
evolve on the ambipolar diffusion timescale tAD ~ 10 i// (Ciolek & 
Basu 2001; tff is the free-fall time), while observed cores have signifi- 
cantly shorter lifetimes ~ 4tff at n H2 ~ 10 4 cm~ 3 (cf. Jessop & Ward- 
Thompson 2000). On the other hand, isolated starless cores appear 
to live £ 3 times longer than typical density fluctuations generated 
by turbulence (cf. Ballesteros-Paredes et al. 2003.) It is possible 
that more elaborate ambipolar diffusion models, incorporating the 
effects of a non-static, turbulent magnetic field in the outer parts 
of the cores and in the ambient cloud, would be more satisfactory 
and could also account for the filamentary shapes often seen on large 
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( > 0.25 pc) scales (cf. Jones & Basu 2002). 

3. VELOCITY STRUCTURE OF CLASS ENVELOPES 

Another approach to constraining the collapse initial conditions 
consists in observing the detailed properties of very young accreting 
protostars (Class objects - e.g. Andre, Ward-Thompson, Barsony 
2000) that still retain detailed memory of their genesis. There have 
been very few quantitative studies of the velocity structure of the 
envelopes of Class objects so far. Here, we describe two recent, 
contrasted examples. 

The first object, IRAM 
04191+1522 (IRAM 04191 for 
short), is relatively isolated. It 
was found at 1.3 mm by Andre, 
Motte, & Bacmann (1999) in 
the southern part of the Tau- 
rus cloud (d = 140 pc). IRAM 
04191 is associated with a promi- 
nent, flattened protostellar enve- 
lope, seen in the (sub)millimeter 
dust continuum and in dense 
gas tracers such as N 2 H + , C3H2, 
H 13 CO+, and DCO+ (Belloche 
et al. 2002). All of the maps 
taken at the IRAM 30m tele- 
scope in small optical depth lines 
show a clear rotational velocity 
gradient across the envelope of 
~ 9kms _1 pc _1 (after deprojec- 
tion), perpendicular to the out- 
flow axis. This gradient is one 
order of magnitude larger than 
those typically observed in star- 
less cores (e.g. Goodman et al. 
1993). Furthermore, the rotation 
of the protostellar envelope does not occur in a rigid-body, but dif- 
ferential, fashion: the inner ~ 3500 AU-radius region rotates signifi- 
cantly faster than the outer parts of the envelope, as indicated by the 
characteristic "S" shape of the position- velocity diagrams obtained 
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Fig. 2. Position-velocity diagram 
along the major axis of the IRAM 
04191 envelope (i.e., perpendicular 
to the flow) based on a C 34 S(2 - 1) 
map taken at the IRAM 30m tele- 
scope (Belloche et al. 2002). Con- 
tours: 0.2 to 0.8 by 0.2 K. The 
dots with error bars mark the ob- 
served velocity centroids. The solid 
curve shows the profile of a model 
with differential rotation [V^(r > 
3500 AU) oc r -1,5 ]. 
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Fig. 3. Rotational velocity (a) and infall velocity (b) inferred in the IRAM 
04191 envelope based on radiative transfer modeling of multi-transition CS 
and C 34 S observations with the IRAM 30m telescope (Belloche et al. 2002). 
The shaded areas show the estimated domains where the models match the 
observations. 



by Belloche et al. (see, e.g., Fig. 2). Fast, differential rotation is ex- 
pected in protostellar envelopes because of conservation of angular 
momentum during dynamical collapse. In the present case, however, 
the dramatic drop in rotational velocity observed at r > 3500 AU 
(Fig. 3a), combined with the flat infall velocity profile (see below), 
points to losses of angular momentum in the outer envelope (see dis- 
cussion in §4 below). 

Direct evidence for infall motions over a large portion of the enve- 
lope of IRAM 04191 is observed in optically thick lines such as CS(2— 1), 
CS(3-2), H 2 CO(2i2 - In), and H 2 CO(3i 2 - 2 n ). These lines are double- 
peaked and skewed to the blue up to ^ 40" from source center, which is 
indicative of infall motions up to a radius Ri n f ^ 5000 AU (cf. Evans 1999, 
2003). Radiative transfer modeling confirms this view, suggesting a flat, 
subsonic infall velocity profile (Vi n f ~ O.lkms -1 ) for 3000 ^, r ^ 11000 AU 
and larger infall velocities scaling as Vi n f oc r -0 5 for r ^ 3000 AU (see Fig. 
3b and Belloche et al. 2002 for details). The mass infall rate is estimated to 
be M in} ~2-3xaf/G~3x 10" 6 M Q yr~ 1 (where a s ~ 0.15 - 0.2 km s" 1 
is the sound speed at T ~ 6 — 10 K), roughly independent of radius. 

Another Class object whose kinematics has been quantified in de- 
tail is IRAS4A in the NGC 1333 protocluster (Di Francesco et al. 2001). 
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Using the IRAM Plateau de Bure interferometer, Di Francesco et al. ob- 
served inverse P-Cygni profiles in H2CO(3i2 — 2n) toward IRAS4A, from 
which they derived a very large mass infall rate of ~ 1.1 x lO -4 M yr -1 
at r ~ 2000 AU. Even if a warmer initial gas temperature (~ 20 K) than in 
IRAM 04191 and some initial level of turbulence are accounted for (see Di 
Francesco et al. 2001), this value of Mi n f corresponds to more than ~ 15 
times the canonical a\jf/G value (where a e ff 0.3 kms -1 is the effective 
sound speed). This very high infall rate results both from a very dense 
envelope (a factor ~ 12 denser than a SIS at 10 K - see Motte & Andre 
2001) and a large, supersonic infall velocity (~ 0.68 km s" 1 at ~ 2000 AU 
- Di Francesco et al. 2001). Evidence for fast rotation in the IRAS4A 
envelope, producing a velocity gradient as high as ~ 40kms~ 1 pc~ 1 , was 
also reported by Di Francesco et al. (2001). 

4. CONCLUSION: COMPARISON WITH COLLAPSE MODELS 

In the case of isolated dense cores such as those of Taurus, the SIS 
model of Shu (1977) approximately describes some global features of 
the collapse (e.g. the mass infall rate within a factor ~ 2) but fails 
to reproduce the observations in detail. The extended infall veloc- 
ity profiles observed in prestellar cores (Lee et al. 2001) and in the 
very young Class object IRAM 04191 (Belloche et al. 2002 - see 
§ 3) are indeed inconsistent with a pure inside-out collapse picture. The 
shape of the density profiles observed in prestellar cores are well fitted by 
purely thermal Bonnor-Ebert sphere models, but the absolute values of the 
densities are suggestive of some additional magnetic support (§2.1). The 
observed infall velocities are also marginally consistent with isothermal col- 
lapse models starting from Bonnor-Ebert spheres (e.g. Foster & Chevalier 
1993, Hennebelle et al. 2003), as such models tend to produce somewhat 
faster velocities. This suggests that the collapse of 'isolated' cores is essen- 
tially spontaneous and somehow moderated by magnetic effects in (mildly) 
magnetized, non-isothermal versions of Bonnor-Ebert cloudlets. Indeed, 
the contrast seen in Fig. 3 between the steeply declining rotation velocity 
profile and the flat infall velocity profile of the IRAM 04191 envelope beyond 
~ 3500 AU is very difficult to account for in the context of non-magnetic 
collapse models. In the presence of magnetic fields, on the other hand, the 
outer envelope can be coupled to, and spun down by, the (large moment of 
inertia of the) ambient cloud (e.g. Basu & Mouschovias 1995). 

Based on a qualitative comparison with the ambipolar diffusion models 
of Basu & Mouschovias, Belloche et al. (2002) propose that the rapidly ro- 
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tating inner envelope of IRAM 04191 corresponds to a magnetically (slightly) 
supercritical core decoupling from an environment still supported by mag- 
netic fields and strongly affected by magnetic braking. A magnetic field 
~ 60 fiG is required at 3500 AU where n H2 ~ 1 — 2 x 10 5 cm -3 , which is 
comparable to the field strengths recently estimated at such densities by 
Crutcher et al. (2003) in three prestellar cores. In this picture, the inner 
~ 3500 AU radius envelope of IRAM 04191 would correspond to the effec- 
tive mass reservoir (~ 0.5 M & ) from which the central star is being built. 
Moreover, comparison of these results with the rotational characteristics of 
other objects in Taurus (Ohashi et al. 1997) indicates that IRAM 04191 
behaves in a typical manner and is simply observed particularly soon after 
protostar formation. The IRAM 04191 example therefore suggests that the 
masses of stars forming in clouds such as Taurus are largely determined by 
magnetic decoupling effects. 



In protoclusters such as NGC 
1333, by contrast, the large over- 
density factors measured in Class 
envelopes compared to hydro- 
static isothermal structures (Motte 
&; Andre 2001), as well as the fast 
supersonic infall velocities and very 
large infall rates observed in some 
cases (§3), are inconsistent with self- 
initiated forms of collapse and re- 
quire a strong external influence. 
This point is supported by the re- 
sults of recent SPH simulations by 
Hennebelle et al. (2003, 2004). 
These simulations follow the evolu- 
tion of a Bonnor-Ebert sphere whose 
collapse has been induced by an 
increase in external pressure P ex t- 
Large overdensity factors (compared 
to a SIS), together with supersonic 
infall velocities, and large infall rates 
( ^ 10a 3 3 /G) are found near point 
mass formation when (and only 
when) the increase in P ext is strong 
and very rapid (e.g. Fig. 4), result- 
ing in a violent compression wave (cf. 




Fig. 4. Density profile (solid curve) 
obtained near point mass formation 
in SPH simulations of the collapse 
of an initially stable, rotating ([3 = 
E rot /E grav = 2%) Bonnor-Ebert 
sphere (T = 10 K) induced by a very 
rapid increase in external pressure 
(with Pext /Pext = 0.03 x the initial 
sound crossing time) (Hennebelle et 
al. 2004). Note the large overden- 
sity factor compared to the p oc r~ 2 
profile of a SIS at 10 K (dotted line). 



Whitworth & Summers 1985). Such 
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a violent collapse in protoclusters may be conducive to the formation of 
both massive stars (through higher accretion rates) and multiple systems 
(when realistic, non-isotropic compressions are considered). Future high- 
resolution studies with the next generation of (sub)millimeter instruments 
(e.g., ALMA - Wootten 2001) will greatly help test this view and shed 
further light on the physics of collapse in cluster-forming regions. 
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